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Regularity of a free boundary with 
application to the Pompeiu problem 

By Luis A. Caffarelli, Lavi Karp, and Henrik Shahgholian * 

Abstract 

In the unit ball i?(0, 1), let u and Q (a domain in W N ) solve the following 
overdetermined problem: 

Au = xn in 5(0,1), G dQ, u = |Vu| = in B(0, 1) \ ft, 

where xsi denotes the characteristic function, and the equation is satisfied in 
the sense of distributions. 

If the complement of ft does not develop cusp singularities at the origin 
then we prove <9ft is analytic in some small neighborhood of the origin. The 
result can be modified to yield for more general divergence form operators. 
As an application of this, then, we obtain the regularity of the boundary of a 
domain without the Pompeiu property, provided its complement has no cusp 
singularities. 

1. Introduction 

In this paper we study the regularity properties of solutions to a certain 
type of free boundary problems, resembling the obstacle problem but with 
no sign assumption, i.e., with no obstacle. Mathematically the problem is 
formulated as follows. Let ft C M. N and suppose there is a function u, solving 
the following overdetermined problem 

(1.1) Au = xn in B, u = |Vu| = in B \ ft, 

where B is the unit ball. 

The question we ask is whether dft is smooth. Indeed, if 9ft is an analytic 
surface then by the Cauchy-Kowalewski theorem, we can always solve the 
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above overdetermined problem locally. We thus ask the reverse of the Cauchy- 
Kowalewski theorem. 

The problem also has a potential theoretic interpretation which is as fol- 
lows. Denote by U the Newtonian potential of f2 (bounded set) with constant 
density (i.e., the convolution of the fundamental solution with xn) and with 
G dfl. Suppose there exists a harmonic function w in B(0,r) (r small) such 
that w = U in Q c (the complement of Q); observe that U is harmonic in Q c . 
This property is called harmonic continuation. 

Next, up to a normalization constant, w — U satisfies equation (1.1). Once 
again the question is whether possession of such a property, harmonic continua- 
tion for S7, will result in the regularity of d£l near the origin. For the interested 
reader we refer to [Sal-3] for similar types of problems. 

It is noteworthy that this kind of problem has been in focus of attention 
in mathematical physics, especially in geophysics [St], [Ma], and in inverse po- 
tential theory [II]. Also, newly developed problems in operator theory reduce 
the analysis of the spectrum of certain hyponormal operators to the study of 
the solutions in problem (1.1), in the complex plane [MY], [P]. We also refer 
to an excellent book by H. S. Shapiro [Shap2], for some basic aspects. 

In order to state our main results, let us define a local solution. 

Definition 1.1. We say a function u belongs to the class P r (z,M) if u 
satisfies (in the sense of distributions): 

(1) An = xn in B r (z), 

(2) u = |Vu| = in B r (z) 

(3) |M|ocvB,(z) < M , 

(4) z £ an. 

We also denote by -Poo(0, M) "global solutions" with quadratic growth, i.e., 
solutions in the entire space M N with quadratic growth \u(x)\ < M(\x\ 2 + 1). 

Remark. If u £ P r (z,M), then 

(1) u{z + x) £ P r (0,M), 

(2) u(z + rx)/r 2 G Pi(0,M/r 2 ), 

(3) Also if ||Aju|| < M, then u(z + rx)/r 2 G Pi(0, M). 

Obviously (1) implies that the class is point independent, so we can always 
consider the class P r (0, M). Our class P differs from that of [Ca2] in two ways. 
We do not restrict the function u to be nonnegative and we replace the uniform 
C 1 ' 1 -norm with a uniform C°-norm. 

This new feature introduces new difficulties, and as a first task we have 
to cope with the optimal regularity of the function u itself before attacking 
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the regularity problem for the boundary dCl. The main tools in our study 
of problem (1.1) will be a monotonicity lemma due to [ACF]; see Lemma 2.1 
below. 

Definition 1.2 (Minimal Diameter). The minimum diameter of a bounded 
set D, denoted MD(D), is the infimum of distances between pairs of parallel 
planes such that D is contained in the strip determined by the planes. We also 
define the density function 

. . , MD({u = | Vu\ = 0} n B(0, r)) 

d r (u) = . 

r 

Now we can state our main results. In Section 7, we apply Theorems I III 
to obtain the regularity of a domain without the Pompeiu property (see §7 for 
details). 

Theorem I. There is a constant C\ = C\(N) such that if u € P\(z,M) 

then 

sup \\Diju\\ < C\M. 

B(*,l/2) 

Theorem II. Let u £ PooiOjM^ and suppose \u — 0} has nonempty 
interior, or S r (u) > for some r > 0. Then u > in and D ee u > in CI, 
for any direction e, i.e. Cl c is convex. Moreover z/liniR—^oo 5r(u) > then u is 
a half -space solution, i.e., u = (max(xi, 0)) 2 /2 in some coordinate system. 

Theorem III. There exists a modulus of continuity a (er(0 + ) = 0) such 
that if u £ -Pi(0, M) and S rQ (u) > cr(ro) for some r$ < 1, then dCl is the graph 
of a C 1 function in B(0, cor^). Here cq is a universal constant, depending only 
on M and dimension. 

Subject to the condition in Theorem III, the analyticity of the free bound- 
ary now follows by classical results [KN], [12]. We thus have the following 
corollary. 

Corollary to Theorem III. Under the thickness condition in Theorem 
III, the free boundary in (1.1) is analytic, in some neighborhood of the origin. 

Theorems I— III are known in 2-space dimensions [Sa2-3]. Indeed, M. Sakai 
[Sa2-3] gives a complete description of the boundaries of all such domains in 
M?. Recently a different approach to this problem was made by B. Gustafsson 
and M. Putinar [GP], where they proved that dCl in (1.1) is contained in an 
analytic arc in IR 2 . This question is still open in higher dimensions. For the 
case u > (the original obstacle problem) the first author has recently proved 
the following: 
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Structure of the singular set ([Ca4]). Let u e Pi(0, M) and u > 0. 
Let y also be a singular point of the free boundary in (1.1); i.e., the free boundary 
does not satisfy the condition in Theorem III, near y. Then there exists a 
unique nonnegative quadratic polynomial (and a unique matrix A y ) 

1 rp 

Qy = 2 X ^y x 

with traceA^ = 1 and such that 

(1) \(u — Q y )(x)\ < \x — y\ 2 a(\x — y\), for some universal modulus of continuity 
a, depending on M only. 

(2) A y is continuous on y. 

(3) // dim(kernel(Aj / )) = k then there exists a k-dimensional C 1 manifold 
Ty, u , such that 

S u n B(y,r) C T y>u , 

for some small r. Here S u indicates the singular points of the free bound- 
ary, i.e., points which do not fall under the hypothesis of Theorem III. 

This fact with no positivity assumption, is studied in a forthcoming paper 
by the first and the third authors. In this paper, however, we will prove 
the analyticity of the free boundary only in the case of "thick" complement 
described in Theorem III. 

Plan of the paper. Section 2 is devoted to some technical tools, which are 
known, but probably not well known in the context used in this paper. In 
Section 3 we carry out the proof for Theorem I, using the monotonicity lemma 
(Lemma 2.1) and the blow-up technique. 

In Section 4 we introduce further lemmas, which will somehow exhaust 
properties of the monotonicity lemma. These are used to prove Theorems II- 
III in Sections 5-6, respectively. In Section 7 we generalize the monotonicity 
lemma to yield for divergence type operators 

Di(aijDju) + a(x)u. 

As a result we obtain the regularity of a domain without the Pompeiu property, 
which we explain to some extent in Section 7. 

2. The monotonicity formula 

In this section we will gather all basic tools used to prove Theorems I— III. 
A fundamental tool, however, is the following monotonicity lemma. 

Lemma 2.1 ([ACF]). Let h\, /i2 be two nonnegative continuous sub- 
solutions of An = in B(x°,R) (R > 0). Assume further that h\h2 = 
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and that hi(x°) = h 2 (x ) = 0. Then the following function is monotone in r 
(0 < r < R) 



(2.1) 



<P(r) = -Af 



|V/ii| 2 \ / f \Vh 2 \ 2 




r) 1% X®\^ 2 I \jB(x°,r) \% X®\^ 



For a fixed direction e, set 

(D e u) + = max(D e u, 0), (D e u)~ = — min(D e u, 0). 

Then we have the following lemma. 

Lemma 2.2. Let u £ Pi(0,M) and consider the monotonicity formula 

(2.1) for the nonnegative subharmonic functions (Deu)^ where e is any fixed 
direction, and denote this by (p(r,D e u). Then 

(2.2) <p\r, D e u) > \{r, D e u)( 7 (T + ) + 7 (r_) - 2), 
where 

(2.3) 7(r±)(7(r±) + iv-2) = A(r ± ), 



with 

(2.4) A(r±) = inf 



/ r± \V'w\ 2 da" 



Here, V is the gradient on S^ -1 , T± = T±(r) is the projection of dB(0,r) n 
{(D e n) ± > 0} onto the unit sphere, and the infimum has been taken over all 
nonzero functions with compact support in F±. 

For a proof of this lemma see (the proof of) Lemma 5.1 in [ACF]. The "set 
function" "/(E) as a function of E C S^ -1 is called the characteristic constant 
of E; see, for example, the paper by Friedland and Hayman [FH]. A result 
of Sperner [Spn] states that, among all sets with given (N — l)-dimensional 
surface area soon on the unit sphere S^ -1 in R N , a spherical cap, i.e. a set of 
the form —1 < c < x\ < 1, has the smallest characteristic constant j(s,N), 
where u>n is the area of the unit sphere in M. N , < s < 1, and c and s are 
coupled by the relations 

s = / (smt) N - 2 dt, c = cos6> , < < it. 

UN Jo 

It thus suffices to consider the function 7 as a function of s = &re&(E)/ujN, 
with s as above. From this [FH, Thm. 2] deduce that for fixed s, 7 is a 
monotone decreasing function in TV (the space dimension), and that the limit 
exists as N tends to infinity 

y(s,N) > 7(5,00). 
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On the other side as N — > oo we will have, for some h > — oo, 

1 f°° 

(2.5) s = ^r ex P (-* 2 /2)*; 

V27T A 

see e.g. [FH, p. 149]. 

Hence for any set £7 £ S w_1 and s as above 

(2.6) 7 (£7, TV) > 7 (£*, TV) = 7 ( a , /V) > lim 7 (s, iV) = 7 (s, oo), 

A^oo 

where E* is the above described symmetrization of E. 
Now from (2.3) (cf. [ACF, p. 441]) one has 

(2.7) T= _l_ +0( jr*). 

In (2.7) if we let N tend to infinity we will obtain, by some tedious calculations 
(cf. [FH, p. 149]), 

(2-8) -^-^ + 0(N~ 3 ) - A, 

where —A is the first Dirichlet eigenvalue of the one-dimensional Ornstein- 
Uhlenbeck operator A — x • V on the set (h, oo), with h as in (2.5). Now 
(2.2)-(2.8) imply that 

(2.9) tp'(r, D e u) > -<p(r, D e u)(A(h+) + A(/i_) - 2), 

r 

where h± are the corresponding constants, in (2.5), for r±. Obviously h + + 
h- > 0. Also, by results of Beckner-Kenig-Pipher [BKP] (cf. also [CK, §2.4]) 
A is convex and A(0) = 1; hence 

(2.10) A(h+) + A(/i_) > 2A(a) = 2. 

In [BKP] it is actually proved that the convex function A satisfies 

A"(0) =4(l-ln2)/7r > 0. 

For convenience we now set 

7(r)=7(r + (r)) + 7 (r_(r))-2. 

Then inserting (2.10) in (2.9) we obtain the following lemma. 

Lemma 2.3. There holds 7 (r) > for all r. Moreover the strict inequality 
holds unless r^_(r) are both half -spheres. In particular if any of the T±(r) 
digresses from being a half -spherical cap by an area-size of e, say, then 

(2.11) 7 (r) > Ce 2 . 
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Remark. The reader may verify by elementary calculus that for E C S^ -1 

with 

area(£) = (- - e)u N , 

i.e., an e digression from the half-spherical cap, we have 

1 1 f°° 1 1 f h 1 

Hence e ~ /i and therefore in Lemma 2.3 we will have 

A(h+) + A(/i_) - 2 > CA"(0)(f£ + /il) « C{h\ + h 2 _) w Ce 2 , 
which gives (2.11). 



3. Proof of Theorem I 

First we need some definitions and notations. 
Definition 3.1. Set 

Sj(,2, m) = sup 

B(«,2-i) 

and define M(z, u) to be the maximal subset of N (natural numbers) satisfying 
the following doubling condition 

(3.1) 4S ij+1) (z,u) > Sj(z,u) for all j eM(z,u). 

Our aim is to prove that Sj < C2~ 2 i , for all j G M(z, u) and some positive 
constant C. An important observation at this point is that if M is empty then 
we may easily (by iteration) obtain the desired estimate. Hence from now on 
we assume 

Lemma 3.2. Let u e Pi(z, M). Then there exists a constant Co = Cq{N) 
such that 

(3.2) Sj(z,u) < C M2- 2j for all j G M{z,u). 

Proof. By the remark following Definition 1.1 we may assume that z is 
the origin. Suppose the conclusion in the lemma fails. Then there exist {uj}, 
{kj} such that 

(3.3) S kj (0, Uj ) > j2~ 2 \ 
with kj e M(0,«j) / 0. 
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Now define Uj as 



. , Uj{2 k ix) „. 
Uj[x) = JV ' in B(0, 1). 



S k . +1 (0,u 3 , 
Then Uj satisfies the following properties: 

2~ 2k i S k .(0, Uj ) 4 

(3.4) WAujW^b < — — < , q j{ '\ < ~ - 0, 

where the second inequality follows from (3.3), and the last inequality follows 
from (3.1), 

(3.5) sup \v,j\ = 1, 

5(1/2) 

(3-6) ll«jl|oo,B < „ ~~F7TT~\ < 4 > 

(3.7) u j (0) = |V« j |(0) = 0. 

Now by (3.4)-(3.7) we will have a subsequence of Uj converging in C 1,a (B) 
(see [GT]) to a nonzero harmonic function no, satisfying uo(0) = |Vuo|(0) = 0. 
For any fixed direction e define 

v = D e uo, Vj = D e Uj, Vj = D e Uj. 

Then, for a subsequence, Vj converges in C 1,a (-B) to v, where Av = 0. Now 
according to Lemma 2.1 (since Vj is harmonic in J7j) 

(3.8) -L f \Vv+\ 2 [ \Vv~\ 2 <C for all r,j, 
r2N Jb(o,v) Jb(o,v) 3 

where C depends on the W 2,2 norm of Uj over the unit ball. By elliptic 
estimates this is uniformly bounded for all j. Making change of the variables 
in (3.8) and letting r = 2~ k j , we will obtain 



L 



\Vv+\ 2 |V^"| 2 < C ( ^ ) <C[^-—\ for all j, 

B(0,1) JB(0,1) K^kj+lJ V z kj J 

where in the last inequality we have used (3.1). Here, and in the sequel, C 
is a generic constant. Next, invoking the Poincare inequality we may reduce 
(improve) the above to 

/ \v+ - M+\ 2 [ \vj - MJ\ 2 < C (^-) for all j, 

JB(0,1) Jb(0,1) \ b kj J 
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where is the mean- value for the functions v^ 1 , on the unit ball. From here 
using (3.3) we obtain, by letting j tend to infinity, 



(3.9) f \v + -M+| 2 f 



\v~ - M~\ 2 = 0, 



/B(0,1) JB(p,l) 

where is the corresponding mean- value for v^. Obviously (3.9) implies 
that either of is constant. Since also v(0) = 0, the constant must be zero. 
In particular v does not change sign. But then the maximum principle gives 
v = 0, i.e., D e UQ = 0. Since e is arbitrary we also have uo is constant. Next 
using uo(0) = we will have uq = which contradicts (3.5). This proves the 
lemma. □ 

Next we will complete the chain in N for the estimate (3.2). 

Lemma 3.3. Let Co = Cq{N) be the constant in Lemma 3.2. Then 

Sj < 4C M2- 2j for all j G N. 

Proof. Let u £ Pi(z, M). Then obviously Si < M. Let j > 1 be the first 
positive integer such that the statement of the lemma does not hold, i.e., 

(3.10) Sj > 4C M2~ 2j . 
Then 

Sj-i < 4C M2- 2 ( J - 1 ) = 16C M2' 2j < 4Sj. 

Hence j - 1 G M(u). By Lemma 3.2, then, 

Sj < Sj-! < C M2- 2( -i-^ = 4C M2- 2 i, 

which contradicts (3.10). The result follows. □ 

From Lemma 3.3 we infer a uniform C 1 ' 1 estimate for the class Pi(z, M). 

Proof of Theorem I. Let u be in Pi(z,M) and set d(x) = dist(x, dft). 
Then by Lemma 3.3 

(3.11) \u(x)\ < CMd{x) 2 for all x G B(0, 1/2). 
Define now 

. . u(x + yd(x)) „ . 

By (3.11) v is bounded on the unit ball, and by its definition it satisfies Av = 1 
in B(0, 1). Hence by elliptic estimates Dijv(0) = Diju(x) is uniformly bounded 
(independent of x). This gives the result. □ 
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4. Further auxiliary lemmas 

In this and the next sections we will frequently use the blow-up of func- 
tions; i.e. for a given u we consider u r (x) = u(rx)/r 2 and let r tend to zero, 
through some subsequence. It is, however, not clear whether the blow-up (the 
limit function) will not be the zero function. Indeed, if u(x) = o(|x| 2 ), then 
any blow-up will be identically zero. To prevent this we need a nondegeneracy 
from below, asserted in the following remark. 

Remark. The function u in (1.1) satisfies 

(4.1) sup u > u(x°) + C N r 2 , 

B(x°,r) 

for all x° G f2. Here CV = l/2iV, if u(x°) > 0, and Cn is somewhat smaller if 
u(x°) < 0. Also r is small enough so that B(x°, r) C B(0, 1). 

The proof of this is given in [Ca2] in the case u > or x° G dQ. The 
general case is proven as follows. Suppose u(x°) < and x° G fl. If there is 
x 1 G B(x , r/2) n dQ, then we apply the above to u in B(x 1 ,r/2) to obtain the 
same estimate with Cn = 1/8N. So suppose the set B(x°,r/2) contains no 
free boundaries, we can apply the mean value theorem for harmonic functions 
to u(x) — \x — x°\ 2 /2N to obtain 

/ (u - u(x )) dx = cr 2+N , 

JB(x°,r/2) 

where c depends on iV only. Prom here one obtains (4.1). 

Now, by (4.1) and Theorem I we may assume that any blow up of functions 
in Pi(0, M) remain in the class. Our next definition will be used frequently 
here and later in Section 6. 

e- close. We say two functions / and g are e-close to each other in a domain 

D if 

sup \f(x) - g{x)\ < e. 
xeD 

Blow-up limit. A blow-up limit uq is a uniform limit on compact subsets 
of R N 

u{rjx) 

uq(x) = hm % — 

where u G P\ (0, M) and rj — > 0. The function u may even change for differ- 
ent j. 

Flat points. We say dil is flat at the origin if there is a blow up uo such 
that the set where uq = is a half-space. 
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Half -space solutions. A half-space solution uq is a global solution that has 
the representation (max(xi, 0)) 2 /2 in some coordinate system. 

For the readers convenience and for future reference we will recall and 
explain some general (known) facts. These will be crucial in the rest of the 
paper. We recommend that a reader unfamiliar with such problems carefully 
verify these facts. 

General Remarks. 

a) By (4.1) and the techniques of [Ca4, Lemma 6] one may show that the 

set 

tt n {x : \Vu(x)\ < e } 

has volume less than Ce, with universal C. One uses only the C 1,1 property 
of the solution, and not the nonnegativity of the function. 
This implies, in particular, that an e-neighborhood 

K £ = {x : dist(x,<90) < e} 

of the free boundary has Lebesgue measure less than Ce, i.e. volume(il" £ ) < Ce, 
with universal C. 

Using covering arguments, we conclude that the free boundary dQ has 
locally finite (N — l)-Hausdorff measure; see [Ca4] for details. 

b) Let Uj be a blow-up of u and suppose Uj converges to uo in C^(M. N ). 
It follows that 

{uo = \Vuq\ = 0} D limjuj = \Vuj\ = 0}, 

(see e.g. ([Ca2], [KS1])); lim denotes the limit set of all sequences {xj},Xj £ 
{ Uj = \ V Uj \ = 0}. 

Next let uo be a blow-up of a sequence Uj and suppose uq = in B(x°, ro) 
for some x°, and ro. Uniform convergence and (4.1) imply that Uj = in 
B(x°,ro/2) for large j. 

A consequence of this is the following 

interior({wo = |Viio| = 0}) C lim{uj = \Vuj\ = 0}. 

c) From a) and b) above we infer an LP convergence of the second deriva- 
tives of Uj to uo; i.e. 

D ik Uj D ik u in Lf oc - norm, 

for 1 < p < oo. This depends on the fact that Auj and Auo differ only 
inside an e-neighborhood of the free boundary d£lo; i.e. on a set of Lebesgue 
measure e. 
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d) From a) and b) we may also deduce that if u G -Pi(0, M) is e-close to 
a half-space solution h = (max(xi, 0)) 2 /2, say, then u = in 

B(0, 1/2) n {xi < -Cy/e}, 

for some constant C > 0. We sketch some details. Let us suppose x° G 
5(0, 1/2) n {zi < 0} n Q. Choose r = |x?|. Then by (4.1) and the closeness of 
u to the half-space solution h we have 

2e > CWr 2 . 



Observe that \u\ < e in 5(x°,r) C {xi < 0}. Hence if x\ < — y / 2e/C l 7v then x° 
cannot be a point of This implies that u = on the set {xi < — y^e/Cjv}. 

e) A consequence of this is the following simple fact: Let u G Pi(0, M) and 
suppose the origin is a flat point with respect to some blow-up sequence u rj . 
Then in B(0, 1/2), u rj is e-close to a half-space solution for small enough rj. 

Lemma 4.1 (essentially due to Spruck [Spk]). Let u G Pi(0, M). T/ien 
any blow-up uq of u is a homogeneous function of degree two, and the set 
{uq = |V«o| = 0} is a cone. 

The proof of Lemma 4.1 in the global case (P<x>(0, M)) is given in details 
in [KS2]. The proof of the local case is similar and therefore omitted. 

Lemma 4.2. Letu e Pi(0, M) and suppose CD\u-u > —e in B(0, 1) for 
some Eq > and C > 0. Then CD\u — u > in B(0, 1/2), provided £q is small 
enough. In particular, if u is close to the half -space solution max(xi,0) 2 /2 in 
B(0, 1), then (by integration and d) in General Remarks) u > in B(0, 1/2). 

Proof. Suppose the conclusion of the lemma fails. Then there is a u G 
Pi(0,M) with 

(4.2) C^m^ ) - u(x°) < 0, 

for some x° G B (0,1/2). Let 

w;(x) = CD\u(x) — u(x) + i^\ x ~ x°\ 2 . 

Then w is harmonic in Q PI B(x°, 1/2), w(a; ) < (by (4.2)) and w > on 
Hence by the maximum principle the negative infimum of w is attained 
on dB(x°, 1/2). We thus obtain 

-e < inf (CDm - u) < — —, 

dB(x°,i/2)nn y ' ~ 87V' 

which is a contradiction as soon as Eq < 1/(8N). The second part, that u > 
in 73(0, 1/2) for it near to half-space solution, follows by d) in General Remarks 
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and integration. The reader should notice that e-closeness for u to a half-space 
solution implies, by elliptic estimates, that the gradient of u also becomes 
(Cy/e) -close to the gradient of the half-space solution. Here C is a universal 
constant. We leave the details to the reader. This proves the lemma. □ 

In the next lemma we will apply the following simple fact, which we for- 
mulate as a remark. 

Remark. Suppose uq is a blow-up solution of u, and <5i(uo) > 0. Then no 
is a degree two homogeneous global solution and the interior of R N \ O is a 
nonvoid cone. 

Indeed by Lemma 4.1, uq has the mentioned properties and R N \ f2 is a 
cone. If the cone has an empty interior, then by General Remark a), M JV \0 has 
Lebesgue measure zero. Therefore by Liouville's theorem uq is a polynomial of 
degree two and hence for all r we will have S r (uo) = 0, which is a contradiction. 

Lemma 4.3. Let u G Pi(0,M), and suppose lim r _^ S r (u) > 0. Then 
u > in some neighborhood of the origin. 

Proof. Let {rj} be a decreasing sequence such that 

lim S r .(u) > 0. 

j— >oc 3 

We blow up the function u through rj to obtain, by Lemma 4.1, a global 
homogeneous solution uq of degree 2 in M N , with 

5i(u ) > 0. 

Hence by the discussion preceding this lemma, ~$l N \ Qq is a nonvoid cone. 

Now choose a direction e and consider the monotonicity formula for D e UQ. 
By degree two homogeneity of uq we have |V-D e uo| is homogeneous of degree 
zero. By scaling, this implies that ip{r, D e uo) must be a positive constant 
unless one of the functions (DeUo)^ 1 is zero. To exclude the first case, observe 
that by Lemmas 2.2-2.3 if ^ then it is strictly monotone since (by the 
condition 5\{uq) > 0) at least one of the sets T± cannot be a half-sphere. In 
the second case we will have D c uq > (or < 0) for any fixed direction e. 

Let x° be a fixed point in and suppose |Vuo(x°)| / 0. Set 

_ Vu (x°) 
V ~ |V« (x°)r 

Then for any directional vector e orthogonal to v we will have D £ uq(x°) = 0. 
Moreover, D c uq > in fio (or < 0). It is harmonic there and it takes a 
local minimum (or maximum). Hence by the minimum (or maximum) prin- 
ciple D e uo = in f^o- This implies that uq is independent of the directions e 
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orthogonal to v and depends only on the direction v. Therefore no is one- 
dimensional in each connected component. Since the only one-dimensional 
solutions are half-space solutions (in each connected components) there must 
be at most two connected components which are half-spaces. Now the assump- 
tion Si(uo) > implies that there must be at most one connected components 
which is a half-space. 

In particular, near the origin, u is close to a half-space solution. Therefore 
the origin is a fiat point. Hence Lemma 4.2 and part e) in General Remarks 
give the result. This proves the lemma. □ 

Our next result will not be used in this paper, however, we include this 
for future references. First we notify the reader of the following obvious fact. 

Remark. Suppose u £ Pi(0, M) and the origin is a point of zero upper 
Lebesgue density for the complement of 0,. Then any blow-up of n at the origin 
is a polynomial of degree two. 

Lemma 4.4. Let u 6 Pi{0, M) satisfy lim r _o 5 r (u) = and suppose for 
some blow-up sequence with limit no, D e uo / 0. Then there is a constant C e 
such that ip(r, D e u) > C e > for all r < 1. 

Proof. Blow up n to obtain a global solution no, in M N with int(Oo) = 
(see the discussion preceding this lemma). Hence no = P, a degree two polyno- 
mial. For any direction e nonparallel to the kernel of P we have D e uo = D e P 
is a nonconstant linear function. Next 



Scaling and letting r tend to zero we will have (see c) in General Remarks) 



where in the last inequality we have used the fact that D e P is a nonconstant 
linear polynomial. Now (4.5) together with the monotonicity formula gives the 



Remark. Lemma 4.4 can be stated in a more accurate way. We can, with 
suitable choice of e, make the constant C e uniformly bounded from below for 
the whole class Pi(0, M). Here is how: Rearrange the coordinate system such 
that the polynomial P in the proof of Lemma 4.4 has the representation 



<p(r) > 4v / \V(D e u)-\ 2 ! \V(D e u) 

J Br J Br- 



(4.5) 




result. 



□ 




i=i 



i=i 
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where m < N . Let now ej be the standard coordinate system and set 

(eH h e m ) 

e = 



jm 
Then 

2 m 
D e P = —=y^aiXi, 



m 

i=i 
and 

m . m 



m ^— ' m z ^— ' TD Z /V z 

Hence for e as above 



i=i i=i 



<p(0, D e u) > 



vol: //, x 2 

2iV 2 



5. Proof of Theorem II 

We remark that in the definition of the global solutions Poo(0, M), we 
require that the functions have quadratic growth with uniform constant M. 
It is noteworthy that this restriction is not superfluous as there are examples 
of solutions to (1.1) in the entire space with M> N \ Q nonvoid and with u of 
polynomial or even exponential growth; see [Shapl]. 

However, we will use Theorem II in connection with blow-up functions, 
i.e., we consider blow-up of functions in Pi(0, M) and these, by Theorem I, 
have quadratic growth near the origin. Hence any blow-up of such functions 
will also be of quadratic growth in the entire space. 

The reader should observe that, by the assumption 5 r (u) > for some 
r > in Theorem II, and the discussion in the remark preceding Lemma 4.2, 
R N \ ft has nonempty interior. It suffices to show that u > 0, since by [Ca2] 
it follows that D ee u > on dQ, and by scaling we also have this property in £1 
(see the details of this scaling argument in [KS2]). We will prove 

Theorem II'. If u £ Poo(0, M) and R N \ has nonempty interior, then 
u>0. 

We split the proof in three cases. 

Case 1. R N \ fl is bounded. 
Case 2. lim^ 00 5^(u) > 0. 
Case 3. lim^oo 5r(u) = 0. 
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Remark. The reader should observe that, by the assumption S r (u) > 
(for some r > 0) in Theorem II, and the discussion in the remark preceding 
Lemma 4.2, M. N \ Vl has nonempty interior. The second part of Theorem II, 
the fact that limJ^(u) > implies u is a half-space solution, is included in the 
proof of Case 2. 

Proof of Case 1. Let U be the Newtonian potential of the complement 
of fl with constant density and such that AU = xn c , i- e - (U = c\.\ 2 ~ N * xn c ', 
for N = 2 we take the logarithmic kernel). Then, since dfl has zero Lebesgue 
measure, A(U + u) = 1 almost everywhere in M. N , and it has quadratic growth. 
Hence by Liouville's theorem it is a second degree polynomial P. 

By translation and rotation we may assume P = ^ a j x2 j + d; observe 
that, by the imposed rigid motion, the origin is not necessarily on dfl any 
more. From here we will have u = P — U in Q. Let now v = x ■ Vu — 2u. Then 
by homogeneity of P — d we will have 

v = —2d — x • VU + 2U -> C asx^oo, 

and v = on dtt. If N = 2 then C = +oo and if N > 3 then C = -2d. 

Suppose first C < 0. Then by the maximum principle v < 0. Now 
by elementary calculus we will imply u{rx)/r 2 is decreasing in r. On the 
unbounded cone- like set K = {rx : x G R N \ r > 1 }, we will have u < 0. 
By subharmonicity, this implies 

= u(x°) < vol(.B) / u, for all x° £ d£l n int(K), 

where -B is a ball in K with center x°. 

Since also u < in K, we will have m = in B, i.e. <9f2 n int(K) = 0. 
Therefore K C M^fl But then R^yQ is unbounded, which is a contradiction. 
Therefore C > 0. 

Now C > implies, by the maximum principle, that u > 0, i.e., u(rx)/r 2 
is increasing in r. Then n < on the truncated cone 

K = {rx: x £ M N \ $7, r < 1}. 

A similar argument as above then shows that u = on K. Hence we conclude 
that the set M. N \ Q has positive Lebesgue density at the origin. 

Now two cases arise: (a) the origin is in the interior of M N \ ft, and (b) 
the origin is on dfl. The first case, along with the monotonicity of u{rx)/r 2 
will result in the positivity of u. As to the second case, we observe that since 
K C {u = 0}, \im r ^o5 r (u) > and by Lemma 4.3, u > near the origin. 
Again by the monotonicity of u{rx)/r 2 in r we will have u > in WL N . □ 
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A simple alternate, but indirect, proof for Case 1 goes as follows: Since 
the Newtonian potential U is a polynomial in O c , it follows from [DF] that Q 
is the exterior of an ellipsoid and from [Shah] it follows that u > 0. 

Proof of Case 2. We first blow up u through a sequence {Rj} (Rj — ► oo) 
such that lim <5r . (u) ^ to obtain a global homogeneous solution uq of degree 
two (see [KS2, Lemma 2.5]). A similar argument as that in the proof of Lemma 
4.3 then implies that no is a half-space solution. In particular this implies 

99(00, D e u) := lim ip(Rj,D e u) = lim ip(l, D e u Rj ) = <p(l, D e u ) = 0, 

Rj— >oo Rj—*oo 

for all directions e. The last equality depends on the fact that a half-space 
solution uq is monotone in every direction, and thus either (D e uo) + = or 
(D e uo)~ = 0. Observe also that we have used the convergence in W 2 ' 2 ; see c) 
in General Remarks. 

Now by the monotonicity formula, 

ip(r, D e u) < <p(oo, D e u) = 0, 

for all vectors e. Hence D e u > (or the reverse); i.e. u is monotone in all 
directions. As in Lemma 4.3 it follows that Vu is parallel at any two points of 
a component of CI, and hence that u is a half-space solution. □ 

Proof of Case 3. Since M. N \ fl is unbounded and lim/j^oo 6r(u) = we 
may conclude that there is a blow-up Uj(x) = u(Rjx)/ R 2 at infinity (Rj — > 
00) with a subsequence converging to a polynomial P in R N , and that P is 
independent of some of the variables. Indeed, by the assumptions in this case 
there is an unbounded sequence x 3 G dQ. Therefore we may take Rj = \x 3 \ and 
obtain that P vanishes, along with its gradient, at the origin and at some other 
point on the unit sphere. Hence by homogeneity (Lemma 4.1) the same is true 
on the whole line generated by these points, and thus the above conclusion. 

Suppose D\P = 0. Then for any point x° £ dil we may consider the 
monotonicity formula (p(r, D\u, x°), which by Lemma 2.1 is nondecreasing in r. 
Hence for r fixed we choose Rj > r to obtain 

<p(r,D lU ,x°) < ip(Rj,D lU ,x°) = ip(l,D lUj ,x°) -► ^(1, £>iP, x°) = 0, 

as Rj — > co. Here we have used that D\Uj — > D\P = 0. This will imply that 
99 = 0, i.e. either of (D\u) + or (D\u)~ is zero. Suppose D\u > (the other 
case is treated similarly). From here we want to deduce that D\\u > 0. To do 
this we set 

— C := inf D\\u, 
which is bounded because u is in Poo(0, M). 



286 



LUIS A. CAFFARELLI, LAVI KARP, AND HENRIK SHAHGHOLIAN 



Let x J be a minimizing sequence, i.e., 

-C = limD u u(x j ). 

j 

We consider a blow up at x J and with dj = dist(x J , dQ). Hence we define 
Uj(*) = - U -j2 " m 5(0,1). 

a 3 

Now by compactness, for a subsequence we will have Uj — > no, where no 
is a global solution. Also DuUj — > DnMo uniformly in S(0, 1/2). The latter 
depends on the fact that Auj = 1 in 5(0,1). Hence we will have a global 
solution no with Ano = Xn where ^o = il(uo). From the minimal properties 
of the sequence x- 7 we also deduce that for all x G B(0, 1/2), 

DhUq{Q) = — C < \\m.DnUj{x) = \\m.Duu(djX + x 3 ) = D\iUq{x). 

j j 

Hence by the maximum principle -Duno = — C in the connected component 
of 0,q that contains the unit ball, we call this 0'. Observe also that the free 
boundary dflo is nonempty. Therefore we assume that z = (z±, • • • , zn) £ Silo- 
Integration gives that in f2 we have the following representation 

D 1 u Q {x) = -Cxi + g(x 2 , • • • , x N ). 

Since also -Di^o = hnij D\Uj > we will have X\ < g(x')/C (of course we 
assume that C > otherwise there is nothing to prove). In particular this 
means that any ray l x emanating at 1 6 SI' and parallel to the xi-axis hits 
<9r2'. Now the component W of 0,q is the under-graph of the function x\ = g/C 
in the xi-direction. Dino > implies that in the negative x\ direction (inwards 
to Q') no is decreasing. Since it is also zero on d£l' it becomes nonpositve in 
W. This contradicts the nondegeneracy (4.1). 
Therefore we conclude C > and hence 

D u u > in Q. 

This along with D\u > implies that n > on lines which hit M. N \ O and 
are parallel to the xi-axis. Our goal will be to prove that for any x in R-^, 

lim u(x\ — m, x 2 , • • • , xn) > 0, 

m—*oo 

which together with D\u > implies n(x) > 0. Now let 

u m (x) = u(xi - m, x 2 , ■ ■ ■ , x N ) (m = 1, 2, • • • ) 

be a family of translations of n. Since the negative xi-axis is in R N \fi (observe 
that D\u > and n(0) = 0) and |n(x)| < M(|x| 2 + 1) we deduce that |n m (x)| < 
CqR 2 on B(0,R) (Cq independent of m). This in particular implies that u m 
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is a bounded sequence; hence there is a converging subsequence with a limit 
function uq in R N . It is then elementary to see that the function uq is also a 
global solution. 

Next using the nonnegativity of D\\u and D\u we infer that the nonneg- 
ative monotone function D\u has a limit at — oo; i.e., 

D\Uq{x) = lim D\u{x\ — m, X2, • • • ,xn) = constant = 0, 

and hence no is cylindrical ({N — l)-dimensional) in the xi-direction. If uq = 
there is nothing to prove; therefore we assume that no is nontrivial and one 
lower-dimensional. Also the set {no = 0} has nonempty interior since it con- 
tains the projection of the interior of the set {n = 0}. But one of the cases may 
occur for this new lower-dimensional function and we repeat the argument (if 
the third case occurs) until we obtain a one-dimensional problem. 

Now the reader may verify, using elementary calculus, that the one-di- 
mensional global solution is nonnegative. This proves n > 0. □ 



6. Proof of Theorem III 

In this section we will give a proof of Theorem III. The proof is based 
on two lemmas that are consequences of already established results, namely 
Theorem II, the techniques in Lemma 4.2, and the main result in [Ca2]. We 
first use Theorem II in conjunction with [Ca2] to obtain uniform flatness for 
the free boundary of u £ -Poo(0, M) provided 8i{u) > e. Next we show that 
for n 6 P±(0,M), there is a uniform neighborhood of such that u > in 
that neighborhood provided n is close to a half-space solution (e close as in 
Section 4). 

The first author's original result for nonnegative solutions in the class 
Pi(0, M) provides us with the following lemma; see [Ca2] or [Ca4]. 

Lemma 6.1. Given a positive number e, there exists t £ such that if 
u S Poo(0, M) and 5\{u) > e, then in B(0,t £ ) the boundary of £1 is the graph 
of a C l function {uniformly for the class) and u is (et £ )-close to a half -space 
solution there. 

Now a simple proof of Lemma 6.1 can be given based on compactness, a 
contradictory argument, and the main result in [Ca2]. 

The proof of the next lemma follows the same lines as that of the proof 
of Lemma 4.2. See also General Remarks. We omit the proof. 
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Lemma 6.2. Let e,s > 0, and suppose u G Pi(0,M) is (es 2 )-close to a 
half -space solution (max(x • e,0)) 2 /2 in B(0,s). Then in B(0,s/2) we have 

(6.1) sD e u-u>0, 

provided e is small enough. In particular, by integration and d) in General 
Remarks, we have u > in B(0,s/4) (again if e is small enough). 

Proof of Theorem III. First we claim that for given e > there exists 

< r £ < t £ , (where t £ is as in Lemma 6.2) 

such that if for u G Pi(0, M) we have S ro (u) > e for some ro < r £ then 
u is (2erQ^)-close to a half-space solution (max(x • e, 0)) 2 /2 in B(0,rot e /2). 
Suppose for the moment that this is true. By Lemma 6.2 it follows that u > 
in P(0, rg/4). Moreover (2erQi 2 ) -closeness to a half-space solution also implies 
(see d) in General Remarks) that u = in 5(0,^/4) n {x\ < — C^/erfy. 
Hence for small e, [Ca2] or [Ca4] applies to conclude that in B(0, coTq), d£l is 
the graph of a C 1 -function with a uniform C 1 -norm. Here cq depends only on 
M and N. 

The modulus of continuity a(r) is then defined by the inverse of the rela- 
tion e — > r £ . 

Now to complete the proof we need to show the (2erQt 2 )-closeness of u to 
a half-space solution in B(0,rot e /2). Suppose this fails. Then for each r,j \ 
there exists Uj G P\(0,M) with 5 rj (uj) > e such that 

sup \uj — h\ > 2er 2 jtl, 

B(0, rj t e /2) 

for all half-space solutions h. 

Set Vj = Uj(rjx)/r 2 -. Then Vj G P 1 / rj (0, M), 5\(vj) > e, and 

sup \vj — h\ > 2et 2 £ , 

B(0,t e /2) 

for all half-space solutions h. 

Now, for a subsequence and in an appropriate space, Vj converges to a 
global solution vo which on one side satisfies the hypotheses of Lemma 6.1, but 
on the other side 

sup l^o — h\ > 2et 2 , 

B(0,te/2) 

for all half-space solutions h. By Lemma 6.1, this is a contradiction. Hence 
the claim holds. 

This completes the proof of the theorem. □ 
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7. Application to the Pompeiu type problem 

Theorems I and III above apply to more general operators of the form 



with dij and a in the C a -class and ay(0) = 5ij (Kronecker's delta function). 
We refer to [Ca3] for these more general operators, and with a(x) = 0. Here 
we treat the case of A + a(x) with a(x) constant. 

To show this we need to make sure that the monotonicity formula works 
for these operators. However, this is not true in general, as one may easily 
observe. On the other hand, the full strength of the monotonicity formula is 
not used in our analysis. 

The monotonicity formula is used either locally near the origin or globally 
when the blow-up functions are considered. For the first case we simply will 
show that the operator L will admit a monotonicity formula which is almost 
increasing. 

The second case is even simpler, as the blow-up of any function in the 
class Pi(0, M, L) (L denotes the dependent on the operator) will converge to 
-Poo(0, M, A), i.e., after rotation the blow-up will satisfy Auo = xn and we 
are back to the same situation as before. 

As to the technique in Lemmas 4.2 (and 6.2) we observe that for the case 
Au + u = 1 we will have A?/; = —CrD e u — 1 + u + cq < if Co is small enough 
and the lemmas work. Spruck's theorem also works with almost no changes. 
We leave the details to the reader. 

A particular problem, where the Laplacian is replaced with the Helmholtz 
operator, is the classical and well-studied Pompeiu problem which can be stated 
as follows. Suppose there exists a bounded domain Q C R N and a function u 
satisfying 



Does it follow that Q is a ball? This question, put forward somewhat differently 
by Dimitrie Pompeiu, has puzzled many mathematicians during the past 50 
years. There are many partial results for this problem, which we will not 
discuss. For background in the Pompeiu problem and for further references 
and how the formulation of Pompeiu is related to the one here see [W] . 

From our perspective it is interesting to know whether the boundary of 
such domains are analytic. Williams [W] proved that if the boundary of £1 in 
(7.1) is Lipschitz then it is analytic. 

It does follow from our results, with the modification below, that under 
the thickness condition of Theorem III, the boundary of such domains are 
analytic. We formulate this in the following theorem. 




(7.1) 



Au + u = xn in R w , 



u = in R \ Q. 
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Theorem IV. Under the thickness assumptions of Theorem III, any do- 
main Q, which admits a solution u to (7.1) has an analytic boundary. 

Now to show the almost monotonicity of the function ip we need to define 
the following functions, 

iv-ui 2 „ s r \vv\ 2 



Then 
(7.2) 



I(r, v) 



ip =- 
r 



I AT— 2 • 



J(r, v) 



B r F 



dB r F 

_ 4+ rI(r,(D e u)+) t rI(r,(D e u)-) 



N-2 - 



J(r,(D e u)+) J(r,(D e u)-) J ' 

Now let F be the fundamental solution of (A + 2). Then 

F = F 1 |x| 2 ~ Ar + F 2 log|:r|, 

where F 2 is zero for odd N, and Fi and F 2 are regular functions not vanishing 
at the origin. Next as in [ACF, Lemma 5.1] (cf. also [Ca3]) we will have, for 
a nonnegative subsolution v to (7.1), 



(7.3) 



J(r, v) < (1 + Cr) [ (A(v 2 /2)F + v 2 F) 

JB{Q,r) 

= (l + Cr) [ (A(v 2 /2)F- V ^AF) 

JB(0,r) z 



= (l + Cr) [ 
Jd 

< (1 + Cr) [ 

Jd 

< (1 + Cr) 



v 



dB(0,r) 



vv nF - —F n 



[v v n r 



2 

2-JV 



dB(0,r) 

rl(r, v) 



+ (N-2) V -r'- N ) 



27(T«(r))' 

where j(T v (r)) is the corresponding r-dependent value for v, as in Lemma 2.2, 
and v n is the normal-directional derivative of v on the sphere dB(0,r). Now 
let us consider the new function 

ip(r, D e u) = ip(r, D e u)e Cr , 

with C to be chosen later so as to make ip nondecreasing. Observe also that 
(Deu)^ 1 are subsolutions to (7.1). Differentiating we obtain 



(7.4) tf(r) 



ip(r) 



rl(r, {D e u)+) + rl(r, {D e u\ 



J(r,(D e u)+) J(r,(D e u)-) 



A + Cr 



Plugging (7.3) in (7.4) we obtain 



tp{r) 



[2 7 (r+)(l + Q{r)) + 2 7 (r_)(l + 0(r)) - 4 + Cr] . 
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Now choosing C large, we see that the extra terms will be taken care of, and 
the result follows as in [ACF, Lemma 5.1]. 
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